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By John R. Spreiter 


SUMMARY 

A method based on assumptions similar to those of Munk's 
airship theory and R.T. Jones* low-aspect— ratio pointed-wing 
theory has been developed to determine simple closed expressions 
for the load distribution, lift, pitching moment, and centeix>f— 
pressure position of inclined slender wing— body configurations 
having flat— plate wings extending along the continuation of the 
horizontal diameters of circular fuselage sections. Expressions 
for the aerodynamic properties of triangular wings in combination 
with conical bodies, semi— inf inite cylindrical bodies, and bodies 
pointed at the nose but cylindrical at the wing root have been 
developed in detail for all ratios of body diameter to- wing span. 

Tn all cases, the lift-curve slope of the wing-body combination was 
less than that of the wing alone . For the case of the triangular 
wing and the body pointed at the nose but cylindrical at the wing 
root, the loss in lift— curve slope reached a maximum of 25 percent 
at the large diame ter-epan ratio of 0.707- With a conical body 
mounted on the same wing, the maximum loss of lift-curve slope was 
only about 8 percent and occurred at about the same diame ter— span 
ratio. 

It is shown that the results are applicable at subsonic and 
transonic speeds, and at supersonic speeds, provided the entire 
wing— body combination lies near the center of the Mach cone. 
Furthermore, it is pointed out that the assumptions related to the 
study of low— aspect— ratio pointed bodies and the study of moderate- 
aspect— ratio pointed bodies traveling at sonic speed both lead from 
Prandtl’s linearized equation for compressible flow to the two- 
dimensional Laplace’s equation in the transverse plane although by 
different means. 


2 


NACA TN No. 1662 


The determination of the potential distribution for an inclined 
moderate— aspect— ratio wing at sonic speed is therefore mathematically 
equivalent to the determination of the potential distribution for an 
inclined low-aspect— ratio wing in an incompressible fluid. 


In the quest for airplane configurations having aerodynamic 
properties favorable for supersonic flight, one of the more 
promising configurations involves the use of a low-aspect— ratio 
wing. When the general layout of such an airplane is considered, 
however, comparatively large fuselages are often found necessary. It 
thus becomes important to study the aerodynamics of a complete wing- 
body combination throughout the entire Mach number range of the air- 
plane. In an incompressible medium, the mutual interference of a 
fuselage and wing of high-aspect ratio (to which lifting-line theory 
is applicable) has been treated by Lennertz, Wieselsberger, Pepper, 
and Multhopp in references 1, 2, 3* and. 4. It is the purpose of this 
note to treat the effect on the aerodynamic loading of the mutual 
interference between a low-aspect— ratio pointed wing and a fuselage 
consisting of a slender body of revolution. 

The aerodynamic properties of slender wing-body configurations 
may be approximated by the method originally used by Munk in studying 
the aerodynamics of slender airships (reference 5) • R. T. Jones 
(reference 6) applied this method to the study of low-aspect— ratio 
pointed wings and Ribner (reference 7) extended it to determine the 
stability derivatives of low-aspectr-ratio triangular wings. The 
essential point in the study of slender bodies by this method is the 
fact that the flow is approximately two-dimensional when viewed in 
planes perpendicular to the direction of motion. Methods of classical 
hydrodynamics may then be employed to determine the load distribution, 
lift, and center of pressure. 

It has been shown by Tsien, Laitone, and R. T. Jones (references 
8, 9, and 6) that the aerodynamic properties of very slender bodies 
of revolution and low-aspect— ratio wings at small angles of attack 
are unaffected by compressibility at subsonic and supersonic speeds. 

A similar result will be shown for slender wing— body combinations. 


INTRODUCTION 
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B 

®b 

®m 

C L 

C La 

C L W 

Cm 

Cm w 

L 

M 

Mo 

S 

U 

w 

X 

a 

c 

c' 

d 

l 

m 


cross-section area of body of revolution (ita 2 ) 
cross-section area of base of body of revolution 


mean cross section of body of revolution 


volume A 
.ength / 


(~ 


lift coefficient l 


lift— curve slope 


dC 

da 


lift coefficient of wing without body 

pitching moment coefficient 
pitching moment coefficient of wing without body 



lift 

pitching moment about apex of wing 
free— stream Mach number 
wing area 
velocity of flight 

complex potential function (cp + i\|r) 
complex variable (y + iz) 
radius of body 
maximum wing chord 

distance from apex to section of maximum span 
semi span of flat plate 
over-ell length of wing-body combination 
additional apparent mass of circular cylinder 
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p static pressure 

q free— stream dynamic pressure 

r ,6 polar coordinates 

s local semispan 

s max maximum semi span 

t time 

v,w velocities in y and z directions 

x, y*z Cartesian coordinates 

x_ ^ distance from apex to center of pressure 

cp velocity potential 

\|r stream function 

a angle of attack 

e downwash angle 

transformed rectangular coordinates 
£ complex variable (r)+iO 

p density of air 

Subscripts 
W wing 

F body 

c compressible 

i incompressible 
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ANALYSIS 

General 

The flow around an inclined wing-body combination of very low 
aspect ratio may "be approx i mated fey considering it to fee two 
dimensional in transverse planes (perpendicular to the fuselage 
center line). It can fee shown as a consequence of this assumption 
that the flow in each transverse plane is independent of that in 
the adjacent planes. Consider a coordinate system moving downward 
through the air with a velocity Ua. The wing-body combination is 
now considered to be flying in the negative x-direction with a 
velocity TJ and angle of attack a so that the fuselage center line 
coincides with the x— axis of the coordinate system and the plane of 
the wing coincides with the z = 0 plane. (See fig. 1(a).) The 
flow pattern, then, in the arbitrary x = x 0 plane during the time 
of the passage of the wing— body combination is approximately similar 
to that of the transverse flow around an infinite cylinder having a 
cross section similar to the local wing-^body section. Observed in 
this plane, the semispan of the wing and the radius of the fuselage 
change with time as the wing— body combination moves through the 
plane. The resulting unsteady nature of the flow pattern produces 
pressure differences between corresponding points on the upper and 
lower surfaces of the wing and fuselag9. The following analysis, 
therefore, consists of three parts: determination of the velocity 

potential for the two-dimensional flow around the wing-body sections, 
determination of the distribution of load over each section, and 
integration of the loading to determine the total lift and pitching 
moment. Several examples are included presenting the total lift, 
center of pressure, and load distribution for typical complete wing- 
body configurations. 


Velocity Potential 

It is necessary for the subsequent analysis to know the velocity 
potential for the unsteady two-dimensional transverse flow field 
around an infinite cylinder, the cross section of which is varying 
with time in such a manner that it always remains similar to the 
wing-body section in the x = x Q plane. Due to the infinite rate of 
pressure propagation in an incompressible fluid, the study of the 
unsteady flow of an incompressible fluid is greatly simplified since 
the flow field at any instant is identical to that of the corresponding 
steady— state flow. The first step in the solution of the present 
problem, therefore, is to determine the velocity potential for the 
steady— state flow around an infinite cylinder having a cross section 
similar to the wing-^body section. In this analysis, only wing-body 
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configurations having circular fuselage sections and flat-plate 
wings extending along the extension of a diameter will be treated. 
The flow around such a section may be obtained from the transverse 
flow around an infinitely long flat plate by application of the 
principles of conformal mapping using the J oukowski transformation. 
Thus we consider the mapping shown in figure 1 in which the £ 
plane will be mapped onto the X plane by the relation 


i 


&2 

= z + x 


( 1 ) 


where 


t = n+H 


and 


X = y+iz 

plane is 

( 2 ) 

plane as opposed 

to the X plane. It is also shown in reference 10 that, if d=2a, the 
flow around a flat plate expressed by equation (2) transforms by 
equation (l) into the vertical flow around a circle of radius a having 
its center at the origin. If d Is taken larger than 2a, the flat- 
plate flow transforms into the desired vertical flow around a cylinder 
consisting of a circular cylinder of radius a with thin flat plates 
extending outward along the extension of the horizontal diameter to a 
distance s from the origin. When the £ plane is transformed into 
the X plane in this manner, the complex potential for the flow in 
the X plane is found to be 


The complex potential function for the flow in the £ 
(see, for instance, reference 10) 

W» = cp* + i\|r' = — iUa V £2-d 2 
where the primed symbols indicate values in the | 


W = 9+ii = — iUa 



2 

-d = 



since the point d in the £ plane corresponds to the point s in the 
X plane. The velocity potential <p for the flow in the X plane may 


then be found by squaring equation (3), substituting X = r (cos 9 + i sin 9), and solving. 
Thus is obtained 



where the sign is positive in the upper half plane (O<0 <n) and negative in the lower 
plane (n<0<2fl). 


Load Distribution 

Once the velocity potential of a flow field is known, the methods of classical hydro- 
dynamics may be applied to determine the pressure at any point in the field. Consider again, 
the case shown in figure 1 where the wing-body combination is piercing the x=Xq plane. As 
previously noted, the flow in the x=Xq plane is considered to be similar to the two- 
dimensional flow surrounding an infinitely long cylinder having the shape of the wing— body 
cross section intersected by the x=x 0 plane. If the radius of the body a and the semispan 
of the wing s are considered to be functions of time, equation (4) may be thought of as 
representing the velocity potential of the unsteady flow in the x=x Q plane . In the case 
of unsteady two-dimensional, potential flow of an incompressible fluid, the pressure at any 
point fixed in the coordinate system is given by (see, for instance, reference 11, p. 19) 

- p = !e + I ( v2+w2 > + *U) (5) 


It may "be seen that thi3 expression reduces to the well-known Bernoulli equation for the 

pressure in a steady flow field when the velocity potential is invariant with time = 0 

ot 

and the arbitrary function of time F(t) is a constant. 
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For any two corresponding points Pi and P 2 (fig. 1(b)), so 
selected that yi = y 2 and z* = — Z 2 , the differential pressure at 
any instant is given by 


_ Pg~Pi 

P P 


8cp2 891 
8t 8t 


i (v 2 2 +W 2 2 ) + i (v X 2 +W! 2 ) 


+ 2 — — (v 2 2 +w 2 2 ) + — (vi 2 +Wi 2 ) 

8t 2 2 


(6) 


&P2 8T1 

since — ^ — = - — by reason of symmetry of the flow field. Now if 
8t 8t 

the points are brought to the wing— body surface 


V 2 2 +W 2 2 = Vi 2 +Wi 2 

the differential pressure between any two corresponding points (or the 
loading) is given by 


Ap 8 <J>l 

— = + 2 — — 
p 8t 


Utilizing the relationship 


(7) 


8cpi cfcpi doc 8cpi 

8t 8x dt 8x 


(8) 


and dividing equation (7) by g 
to be N 

k ^ 

_ Pz~Pu = dx _ 

q " q U 


U 2 , the loading coefficient is found 


/ 8cpi ds 8q>x da\ 

U ^ 8s dx 8a d x J 


(9) 


The load distribution may now be obtained by substituting the 
expression for the velocity potential given in equation (4) into 
equation (9) and letting 6=0 or 0= it for the wing loading and 
r=a for the fuselage loading. The loading over the wing is then 
found to be given by 
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and that over the fuselage is given by 



In Cartesian coordinates, the loading over the fuselage is 


( 11 ) 



Total Lift and Moment 

The total lift and pitching moment of a complete wing-^body 
combination may be determined by integrating the loading over the 
entire plan— form area. It is convenient to carry out the integra- 
tion by first evaluating the lift on one spanwise strip and then 
integrating these elemental lift forces over the length of the wing- 
body combination. The lift on a spanwise strip of width dx is 
given by 


dL •= 2q dx 



<3 /F 


a sin 9 dd + 


i 



(13) 


or, in Cartesian coordinates 






# 


dL = 2q dx 


(14) 
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When the indicated operations are performed, the following expressions 
for the elemental lift on the wing and body are obtained. 



equations (15a) and (15b) may be combined and simplified to give the 
following expression for the total lift on an elemental spanwise strip. 
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The lift, pitching moment, and center of pressure of the complete 
wing-body combination may now be determined by integration of the 
lift of all the elemental strips 



where the integration interval extends from the most forward point 
to the most rearward point of the wing-body configuration. The lift 
coefficient, moment coefficient, and center of pressure may be 
determined from equations ( 18 ), (19), and (20) by division by 
appropriate constants 


C 


L 



( 21 ) 


C 


m “ 



*0 .p . ^m 



( 22 ) 

( 23 ) 


where S is the reference area and c the reference chord or 
length. 


12 


MCA TN No. 1662 


Effect of Compressibility 

In contrast to the well-known infinite-aspect-ratio case where 
the pressures on the surface of a wing are influenced by compressi- 
bility in a manner described by the Prandtl-Glauert relation, it has 
been shown by several investigators that the pressures on very low- 
aspect-ratio wings and very slender bodies of revolution are 
unaffected by compressibility. This result has been found by Jones 
(reference 6) for low-aspect-ratio pointed wings at both subsonic 
and supersonic speeds. B. Gothert (reference 12) extended the low- 
aspect-ratio rectangular wing theory of Bollay (references 13 and 14) 
to include the influence of compressibility and found no effect in 
the subsonic range. For a very slender inclined body of revolution 
a ^ S Q^ S ° niC 311(1 su P ersonic speeds* Laiton© and Tsien (references 9 
and «) have found that the loading was unaffected by compressibility. 
That such is also the case for slender inclined, pointed wing— body 
combinations follows from consideration of the basic differential 
equation of linearized compressible flow. In addition, it will be 
shown that the aspect-ratio range to which the theory is applicable 
becomes larger as the Mach number approaches one. 

Prandtl (reference 15) has found the linearized differential 
equation for the velocity potential of compressible flow to be 


(1-Mq 2 ) |-2 

OX 2 


dfcp 

c>y 2 



(24) 


In the development of the expressions for the forces on long slender 

wing-^body combinations, it has been assumed that — is so much 

ox 2 

, d 2 cp c) 2 © 

smaller than — and — — that the first term of equation (24) 
oy 2 oz 2 


may be neglected. Therefore, so long as the term (l-Mn 2 ) ^ — 

ox 2 

in the differential equation remains small, Mach number will have 
little influence on the distribution of the velocity potential. 
Consequently, Mach number has little effect on the aerodynamic 
characteristics of a long slender wing— body combination at either 
subsonic or supersonic speeds. It is Immediately apparent that the 
Mach number cannot be increased indefinitely, for then the coeffi- 

d 2 cp 

cient of — - becomes so large that the first term will no longer 
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be negligible. The required condition will be satisfied, however, 
if the body has a pointed nose, the wing a pointed plan form, and 
the entire wing-body combination lies near the center of the Mach 
cone. All these conditions, however, correspond to those originally 
assumed in the derivation of the expression for the velocity potential 
(equation (4)). Therefore, the present theory is applicable at 
supersonic speeds, as well as subsonic speeds, provided the entire 
wing-body combination lies near the center of the Mach cone. 


It has been shown by Robinson and Young (reference 16) that, 
for finite aspect ratio, the linearized theory of compressible flow 
(equation (24)) remains theoretically consistent and yields finite 
and continuous lift— curve slopes in the transonic range. Recent 
experiments on triangular wings at transonic speeds support this 
contention by indicating agreement between measured and computed 
lift-curve slopes. Therefore, to predict the flow around a body 
traveling at or very near sonic velocity, it is correct, unless the 


term 


d 2 cp 

- — becomes extremely large, to let 
ox 2 


Mq=1 


and solve the 


remaining equation for the potential distribution. The remaining 
equation is the two-dimensional Laplace's equation in the transverse 
plane. This means that, although the velocity potential may vary 
in the longitudinal direction, its value at each point may be 
determined solely by studying the flow in the transverse plane 
containing the point in question. Therefore, since this is precisely 
the manner in which the potential distribution was obtained, the 
results of the present analysis are applicable at transonic speeds. 

In fact, the present theory is most applicable to wing-body 
combinations of moderate aspect ratio if the Mach number is one, since 

it is then no longer necessary to assume that - — - is very much 

, , , <3 2 cp <3 2 cp ^ x2 

smaller than — - and — - . 

6y 2 dz 2 


In retrospect, the assumptions related to the study of low-aspect— 
ratio pointed bodies and the study of moderate-aspect— ratio pointed 
bodies traveling at sonic speed both lead from the Prandtl equa- 
tion (equation (24)) to the two-dimensional Laplace's equation in 
the transverse plane although by different means . The low-aspect— 


ratio theory neglects the term (l— Mr, 2 ) ^1? 

ox 2 


in comparison with 


a 2 ? 

dy 2 


and 


a 2 ? 

Sz 2 


d 2 9 

because — — is very small; while the moderate-aspect— 
ox^ 


ratio sonic theory neglects the same term because (1 -Mq 2 ) is zero. 
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Hiua the determination of the potential distribution for an inclined 
moderate— aspect— ratio wing at the speed of sound is mathematically 
equivalent to the determination of the potential distrihution of an 
inclined low-aspect-ratio wing in an incompressible fluid. 


EXAMPLES 

For a given wing-tody configuration complying with the general 
requirements of the present theory, the load distribution may be 
determined directly by substituting the proper values for the body 
radius and wing semispan and their rate of change with x into 
equations (10) and (ll) . In addition, closed expressions for the 
lift, pitching moment, and cent© r-of— pressure position of several 
elementary configurations may readily be found by simple integration 
of the integrals indicated by equations (21), (22), and (23). 

Several such examples will be presented in detail in this section, 
and the results will be compared in the following section with those 
obtained from linear theory and from experiment. 


Pointed Low^Aspect-Ratio Wing 

Although the assumptions of this note have been used previously 
by R. T. Jones in reference 6 to determine the aerodynamic properties 
of law-aspect-ratio wings, the load distribution, lift, and 
pitching moment will be rederived for completeness of presentation 
and to show a simple application of the preceding expressions. The 
aerodynamic properties of a low— aspect— ratio wing without fuselage 
may be determined by letting 


By substitution of these values into equation (10), it follows that 
the load distribution along any elemental spanwise strip is 



(25) 


The loading (fig. 2(a)) thus shows an infinite peak along the 
leading edge of the wing. The total load on an elemental spanwise 
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strip is found from equation (17) to be 



, ds 

4rtoqa — 
dx 


( 26 ) 


Equations (25) and (26) show that the development of lift by the 
long slender wing depends on an expansion of the sections in a 
downstream direction. Accordingly, a part of the wing having 
parallel sides would develop no lift, while a part having contract- 
ing width would have negative lift with infinite negative loads 
along l^e edges. In the actual flow, however, as B. T. Jones points 
out preference 6), toe portion of toe wing behind the maximum cross 
section will lie in the viscous or turbulent wake formed over the 
surface ahead. Consequently, toe infinite negative loads will not 
be developed on these edges. With toe aid of the Kutta condition, 
Jones then concludes that no lift is developed on sections aft of 
toe maximum cross section. This is known to be an oversimplifica- 
tion of toe truth and considerable caution should be exercised in 
applying toe present results in toe case of constant or gradually 
contracting width. 


■® ie lifl coefficient for this wing is found by integration of 
toe load on toe elemental strips between toe leading edge and the 
widest section as indicated by substituting equation (26) into 
equation (21) 




^TTLft.Y 


sds = - a 


= - Aa (27) 


where c' is toe effective wing chord and 


S 

dCj, 


= A, toe aspect 


It is seen that the 1 if t — cur v e slope — depends only 

on toe aspect ratio. It should be noted, however, that the actual 
lift force depends only on toe span and angle of attack and not on 
toe aspect ratio or toe area. 

By similar substitution and integration by parts of equation 
(22), toe pitching moment about toe leading edge is 


c * 

n, 1 P , ds n c* 

-m = - — / 4nas — xdx = - - — a 

9cJ dx 2 c 


^• s max 2 


Ms 2 ) 


ml * c' 

( 28 ) 


1 - 


(s 2 ) 

8 max 2 


s] 


1 6 
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1 n ° 1 

where (a 2 ) - — r / s 2 dx and where moments tending to produce 

m 0 Jo 

a nosing-up rotation are considered positive. The center— of— 
pressure location is then found by dividing the moment coefficient 
by the lift coefficient as indicated in equation (23). 


Xc.p. 

Cm c * 

Ms a )m' 

c* 

(s 2 )m~ 

c 

H 

I 

t?l 

II 

°l 

1 SA 


^ax 2 - 


(29) 


For a more specific example, consider a triangular wing moving 

point foremost. Then since (s 2 ) m = and c'=c, the pitching- 

moment coefficient and center-of -pressure position are given, 

respectively, by (^ = - ^ Aa and = | . The center of pressure 

is seen to be at the two— thirds chord point or the center of area. 


Pointed Slender Body of Bevolution 


The present method for treating the flow around long slender 
bodies was introduced by Munk in reference 5 for the determination 
of the distribution of forces along the longitudinal axis of a body 
of revolution (airship hull) . In the present section, these results 
will be rederived. In addition, expressions for the total lift, 
pitching moment, and load distribution will also be presented. 

For the slender pointed body of revolution, the following re la 
tions exist: 

a da _ ds 

s dx dx 

where is not necessarily constant. If these values are 

substituted into equations (ll) and (12), the loading distribution 
along any elemental strip is 

(t) f = 8a- al sln 9 = ^ £> _ h 


(30) 
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The load distribution (fig. 2(b)) is thus seen to be elliptical, 
being zero at the extremities of a horizontal diameter and a maximum 
at the midpoint. The total load on an elemental apanwiae strip is 
found from equation (17) to be 


dx 



= 4noqa — = 2aq — 
dx dx 


(3D 


where B is the local cross-section area. It is seen that equation 
(31) is identical to equation (26) for the integrated load on an 
elemental spanwise strip of a triangular wing, even though the 
distribution of load in the two cases is widely different. In 
contrast, however, to equation ( 26 ), which is to be applied only to 
wings of increasing span, equation (31) may be applied to bodies of 
revolution in regions of either increasing or decreasing radius, 
since the Kutta condition does not apply to bodies of revolution. 
Thus, in general, the lift and pitching moment of a body of revolu- 
tion are different from those of a wing of identical plan form; 
however, if the maximum diameter of the body of revolution is at 
the base station, its lift and pitching moment are equal to those 
of a wing of identical plan form at the same angle of attack. 


As before, the lift coefficient will be determined by substituting 
equation (31) into equation (21) . Taking the area of the base cross 
section B^ as the reference area and integrating over the length 
of the body l the lift coefficient is found to be 



since the cross-section area B is B^ at x=l and zero at x=0. 

It is thus seen that the lift of a slender body of revolution depends 
only on the cross-section area of the base, and is independent of 
the general shape of the body. A possible effect of viscosity is 
indicated by such a relationship since the effective base area of 
the body will be larger than the true base area by an amount 
dependent on the boundary— layer thickness. Therefore equation (32) 
will probably tend to underestimate the true lift— curve slope, 
particularly at lower Reynolds numbers where the boundary— layer 
thickness is greatest. 

By similar substitution and integration by parts, the moment 
coefficient about the leading edge is 
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Cjn - 


-1 

r l dB / 

i-i' 


/ 2a — xdx = -2a ( 

B \>l 

'0 dx ' 


mean 

cross-section area 

(i.e.. 


where % is t 
body divided by the length). The center— of -pressure 
then found through use of equation (23) to be 


(33) 


volume of the 
location is 




( 3*0 


For a more specific example, consider a cone moving point foremost. 
The base cross-section area is 


= «a 2 


The mean cross-section area is 

Bm = j na2 

The center of pressure is thus seen to be at the two— thirds point as 
would be anticipated by the conical nature of the load distribution 
for this case. 


Triangular Wing With Conical Body 


The first example of a wing— body combination to be considered 
is that of a conical body mounted on a triangular wing so that their 
vertices coincide. The geometry of such a configuration requires 
that 

a da/dx 
b ~ ds/dx 


where both — and -r— are constants. If these values are 
dx dx 

substituted into equations (10) and (ll) as described in the two 
preceding examples, the load distribution along any elemental strip 
on the wing is given by 



1+k 4 -2k 4 


s2 

y2 


1+k 



for a <y < s 


(35a) 
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and. on the "body hy 




l+k 2 ) 2 



for o <y <a 


(35b) 


Figure 2(c) shows the load distribution on a typical wing— body 
combination of this type together with the load distribution on the 
same wing without body. 

The integrated load on an elemental strip is 


dx \ a / q 


where 


^ naqs S | 1+k4+ ^ 2k ( 1— k 2 ) — ( 1+k 2 ) 2 s in 1 ^p:] j= 4naqs^(l+E) 

(36) 


E = k4 + 27 


2k(l— k 2 )— (l+k 2 ) sin 1 


-i 2k 


l+k2 


The lift coefficient for the entire conical wing-^body combination is 
then 


C L = | Act, (1+B) = C Lw (1+E) 


(37) 


where Cl^ is the lift coefficient of the basic triangular wing. 


The area and aspect ratio of the wing— body configuration are 
considered to be equal to those of the basic wing. Due to the 
radial nature of the lines of constant pressure, the center of 
pressure lies at the two— thirds chord point 


c 3 

The moment coefficient is then obviously 


(38) 


c m = — | Aa(l+B) = Cj^d+E) (39) 


where, similar to before, C mw . represents the pitching moment of 
the basic wing. Figure 3 shows the variation of C 


i« /c S, 


with 


ratio of body diameter to wing span for this type of wing-body 
combination. While the wing alone and body alone have identical 
lift— curve slopes since the widest section is at the trailing-edge 
the lift— curve slope of the wing-body combination is always less 
than that of either a wing or body alone. The maximum loss of 
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lift— curve slop© (about 8 percent) occurs when the body-radius 
wing— eemispan ratio is approximately 0.7. 


Triangular Wing on a Semi-Infinite Cylindrical Body 

The next example to be considered is that of a triangular wing 
mounted on a semi-infinite cylindrical body. The essential relation- 
ships to be used are that 

da 

— =0 
dx 

and that ds/dx is constant. By using these relationships as in the 
previous examples, it is found that no lift is carried on the body 
ahead of the leading edge of the root chord. Behind this point, 
however, lift is carried on both the wing and body and is distributed 
on any elemental strip of the wing in a manner described by 



and on the body by 



for 


a< y< s 


o< y <a 


(40a) 


(40b) 


The load distribution at one longitudinal station of a typical wing- 
body configuration of the type considered in this example is shown 
in figure 2(d) . For purposes of comparison, the load distribution 
over the same wing without the body is also indicated in figure 2(d) . 


The integrated load on an elemental strip is given by 


d 

dx 



ds 

4KCUJS ^ 



(41) 


By integration along the length of the body, the lift coefficient for 
the complete wing— body combination, based on the area of the basic 
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triangular wing without fuselage is found to be 


C 


L " 


jt 

2 Act 



^ax 2 





It may be seen from equation (42) and figure 3 that the addition of 
a semi— inf inite cylindrical body to a triangular wing produces a 
loss in lift— curve slope just as in the preceding example with the 
conical body. With the cylindrical body, however, the lift-curve 
slope has no minimum value, but continues to decrease as the radius— 
semispan ratio increases until finally, when the latter ratio is one 
(corresponding to a body without wings), the lift-curve slope is zero. 
This is as it should be, since a semi— inf inite cylindrical body has 
zero lift— curve slope. The moment coefficient about the vertex of the 
basic triangular wing is 



(43) 


The center— of -pressure position of the complete wing— body combination 
is given by 



Since the center of pressure of the wing alone is at the two— thirds 
chord point, it may readily be seen the second term of equation (44) 
represents the change due to the addition of the body. Figure 4 
shows the variation of the center— of— pressure position with the ratio 
of body radius to wing semispan. In contrast to the constant center^ 
of— pressure position of the previous example for the triangular— wing, 
conical— body combination, the center of pressure of the triangular 
wing, semi— inf inite cylindrical body combination moves rearward as 
the body radius becomes larger with respect to the wing semispan. 


Triangular Wing on a Pointed Body 

The case of a triangular wing mounted on a pointed body, closed 
in an arbitrary manner at the nose but cylindrical along the wing 
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root chord, may he studied by combining the results of two previous 
examples. The portion of the wing— body combination ahead of the 
leading edge of the wing root may be considered to be equivalent to 
the arbitrary body of revolution treated in the second example. The 
portion of the wing— body combination aft of the leading edge of the 
wing root is equivalent to a triangular wing mounted on a semi- 
infinite cylinder as discussed in the preceding example . The load 
distribution and the integrated load on any elemental spanwise strip 
are then the same as those given in the corresponding example . 

The lift coefficient is found by adding the lift forces of the 
component parts of the wing— body combination and dividing by the 
dynamic pressure q and the characteristic area, again taken to be 
the area of the basic triangular wing. The lift coefficient is then 
found to be 


C L “ § Aa 



a 


3 max* 


+ 

9 max 4 / 



(^ 5 ) 


Figure 3 shows the variation of the lift-curve slope with tody-radius 
wing— semi span ratio, A comparison of the lift— curve slopes shows 
that the loss in the lift of a triangular wing resulting from the 
addition of a tody having a pointed nose is much less than that 
resulting from the addition of a semi— inf inite tody. 

The moment coefficient for this wing— body combination may be 
found in a manner similar to that used in finding the lift coeffi- 
cient, taking care to transfer the moments of both component parts 
to the same axis. The moment coefficient about the vertex of the 
basic triangular wing is 


Cm 


TtctA / . a 3 a 4 \ 

~3~ \ ~ a max 3+ 'hoax 4 ' 


2ita a 3 B m /Z + a A 
S _ s max n v 3 8 max^ 



(l \ a3 , ’ a4 N 

, 2na 

[a 3 B m ( 

(L 
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1 + 1 
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_ 3 max * 


s max>' 
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(46) 


where a represents the radius of the cylindrical portion of the 
fuselage, l the over-all length of the wing-body combination, and 
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Bjjj the mean cross-sectional area (i.e., volume divided by length) 
of the portion of the body ahead of the leading edge of the wing 
root . 


COMPARISON WITH OTHER RESULTS 

As shown in the preceding sections, it is a comparatively simple 
matter to calculate the load distribution, lift, and center of 
pressure of complete wing— body configurations by means of the present 
theory. It has been shown that the theory is most applicable at Mach 
numbers near one or for configurations having very low-aspect— ratio 
wings. Its accuracy at other Mach numbers or at larger aspect ratios 
can best be assessed by comparison with experiment or more nearly 
exact theory, where available. 

Comparisons with available theoretical and experimental lift- 
curve slopes of triangular wings of varying aspect ratio at superb- 
sonic and subsonic speeds are shown in figures 5(a) and 5(b) j 
respectively. In the supersonic range (fig. 5(a)), the linear theory 
solution of Stewart, Brown, and others (references 17 and 18) for 
the variation of lift-curve slope with aspect ratio is shown for 
Mach numbers of 1.0, 1.2, and 1.4. At a Mach number of 1.0, it is 
seen that the present theory exactly predicts the linear theory 
value of the lift-curve slopes of triangular wings of any aspect 
ratio. Increasing the Mach number decreases the degree of correla- 
tion at the larger aspect ratios. In summary, this figure indicates 
that the present theory is very accurate for slender wings at low 
supersonic speeds where the wing is near the center of the Mach 
cone, and decreases in accuracy as the wing becomes larger with 
respect to the Mach cone. 

In the subsonic case (fig. 5(b)), no lifting— surface theory 
for the triangular wing comparable to the supersonic triangular- 
wing theory exists, and all comparisons will be made directly with 
experiment. Three test points from reference 19 are shown for 
wings of aspect ratio 0.5 j 1.0, and 2.0 tested at very low Mach and 
Reynolds numbers in the Langley free— flight tunnel. As in the 
supersonic case, the accuracy is best at very low aspect ratios 
and decreases as the aspect ratio increases. 

A comparison between lift— curve slopes for a complete wing- 
body combination consisting of a conical body and a triangular 
wing calculated by the present theory and by supersonic conical- 
flow theory is shown in figure 3» A curve presented by Browne, 
Friedman, and Hodes (reference 20) for the lift— curve slope of a 
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wing— body configuration consisting of a conical body having a fixed 
radius of 0.1322 the Mach cone radius and a triangular wing of 
varying span is shown by the dotted line in figure 3 together with 
the corresponding curve obtained by the present theory. These two 
curves never differ by as much as 1 percent, indicating that the 
present theory and the conical— flow theory are in close agreement 
in predicting the lift— curve slope at supersonic speeds of a wing- 
body combination consisting of a slender conical body and a low- 
aspect— ratio triangular wing. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 


APPENDIX 

Correction to Loading on Portion of Fuselage 
Aft of Wing Trailing Edge 

A method for the calculation of the aerodynamic loading on the 
entire surface of a slender pointed wing-body combination has been 
presented based on the assumption that the flow in each transverse 
plane is independent of that in the adjacent planes. It was noted 
that the results so obtained were not applicable to the portion of 
a wing situated behind the widest section because the flow in this 
region was influenced to a prohibitive degree by the downwash field 
of the sections further forward. For the same reason, the results 
are also inapplicable to the portion of the fuselage aft of the 
wing trailing edge, particularly when the fuselage diameter is 
small in comparison with the wing span. Since the fuselage is 
usually extended behind the trailing edge of the wing, it is desir- 
able to determine a correction to apply to the loading expressions. 

With assumptions more restrictive than those of the main body 
of this note, it is possible to obtain an estimate of the corrected 
loading on the fuselage afterbody. The necessary assumptions are 
that the downwash velocity in the vicinity of the fuselage with 
fuselage removed is known, and that the downwash velocity remains 
constant throughout the entire transverse plane at each longitudinal 
station. It is immediately apparent that the latter assumption is 
not entirely correct, but it is true that the downwash velocity is 
approximately constant over a region of limited lateral extent at 


NACA TN No. 1662 


25 


each longitudinal station. Since the forces on a body are produced 
predominately by the flow field near the body, this assumption 
should be a valid one as long as the fuselage sections remain in 
regions of relatively constant downwash velocity in each transverse 
plane. This means that the fuselage diameter must be small in 
comparison with the wing span. 


With the foregoing assumptions, the loading on the fuselage 
afterbody may be determined by an extension of the present method. 
Consider, as in figure 1(a), the flow in the transverse plane as 
the fuselage afterbody is piercing the x=xo plane. The flow field 
corresponding to that of figure l(b) would then be that of the 
vertical flow around a circular cylinder. As in the previous 
analysis, the fuselage radius would, in general, appear to be 
varying with time. In addition, since the downwash velocity varies 
with distance behind the wing, the velocity of the vertical flow 
would also appear to be varying with time. The correct expression 
for the lift on each strip across the fuselage may then be obtained 
by substituting the local angle of attack a-€ for the airplane 
angle of attack a in equation ( 31 ) and adding a correction term 
for the effect of the longitudinal gradient of the downwash velocity. 
The latter correction term may be determined very simply using the 
additional apparent mass concept. The correction to the lift force 
on an elemental strip of unit width across the fuselage afterbody 
is then given by 



m — = -2Bq 
dt 


d£ 

dx 


(Al) 


In this equation, the additional apparent mass of a unit length of 
a circular cylinder of cross-section area B is (see, for instance 
reference 11, p. 77) 


m = pB 


(A2) 


and the vertical velocity in any transverse plane is 


W = U ( Ob — € ) 


(A3) 


where e is the downwash angle. The total lift on each elemental 
strip of unit width of the fuselage afterbody is then 



2q 


(ar-e) 


dB de 7 

B — 

dx dx J 


(A4) 


d 

dx 
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Since the downwash angle numerically equals the angle of attack on 
the wing surface immediately ahead of the wing trailing edge and 
decreases in value as the distance from the trailing edge increases, 

it is apparent that € has a positive value and that ^ has a 

cue 

negative value. Thus, effects of downwash angle in equation (A4) 
tend to cancel each other. Another item which should he mentioned 
is that the downwash behind wings varies considerably with Mach 
number. Consequently, compressibility will affect the lift on 
fuselage afterbodies. 

At subsonic speeds, an upwash exists over the portion of the 
body extending ahead of the wing, although this upwash is of 
considerably smaller magnitude than the downwash behind the wing 
for the slender pointed wings considered here. Equation (A4) may 
be applied to determine the magnitude of the corrected loading 
taking into account the upwash. 

By using methods similar to those developed in the main body 
of this note rather than the shorter additional apparent mass 
methods, it can be shown that the load distribution across each 
strip of the fuselage afterbody is elliptic . 
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Figure /.- View of wing- body combination showing coordinate axes and two-dimensional flow 

fields. 
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Loading over wing without body 


Figure 2 . . - Load distributions over spanwise sections of four 
configurations. 


MCA TN No. 1662 


31 



Wing semi span s 


Figure 3.- Lift-curve slope ratios for three wing-body configu- 
rations, 
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Figure 4- Center-of-pressure position of two wing-body configu- 
rations. 
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Figure 5.- Comparison of theoretical and experimental lift-curve slopes for triangular wings. 
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